Astronomy & Astrophysics manuscript no. Thurgood'McLaughlin'2012a 


©ESO 2012 


August 30, 2012 





(N 

o 

< 
ov 

(N 



Or 

6 

o3 



Linear and nonlinear MHD mode coupling of the fast 
magnetoacoustic wave about a 3D magnetic null point 

J. O. Thurgood and J. A. McLaughlin 

School of Computing, Engineering & Information Sciences, Norfhumbria University, Newcastle Upon Tyne, NE1 8ST, UK 
Received 20 June 2012 / Accepted 29 July 2012 

ABSTRACT 

Context. Coronal magnetic null points have been implicated as possible locations for localised heating events in 2D models. We 
investigate this possibility about fully 3D null points. 

Aims. We investigate the nature of the fast magnetoacoustic wave about a fully 3D magnetic null point, with a specific interest in its 
propagation, and we look for evidence of MHD mode coupling and/or conversion to the Alfven mode. 

Methods. A special fieldline and flux-based coordinate system is constructed to permit the introduction of a pure fast magnetoacoustic 
wave in the vicinity of proper and improper 3D null points. We consider the ideal, /? = 0, MHD equations which are solved using 
the LARE3D numerical code. The constituent modes of the resulting wave are isolated and identified using the special coordinate 
system. Numerical results are supported by analytical work derived from perturbation theory and a linear implementation of the WKB 
method. 

Results. An initially pure fast wave is found to be permanently decoupled from the Alfven mode both linearly and nonlinearly for 
both proper and improper 3D null points. The pure fast mode also generates and sustains a nonlinear disturbance aligned along the 
equilibrium magnetic field. The resulting pure fast magnetoacoustic pulse has transient behaviour which is found to be governed by 
the (equilibrium) Alfven-speed profile, and a refraction effect focuses all the wave energy towards the null point. 
Conclusions. Thus, the main results from previous 2D work do indeed carry over to the fully 3D magnetic null points and so we 
conclude that 3D null points are locations for preferential heating in the corona by 3D fast magnetoacoustic waves. 

Key words. Magnetohydrodynamics (MHD) - Waves - magnetic fields - Sun: corona - Sun: magnetic topology - Sun: oscillations 
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1. Introduction 

An abundance of observational data from solar instruments in- 
cluding SOHO (e.g. Ofman et al. 1997; DeForest & Gurman 
1998), TRACE (e.g. Nakariakov et al. 1999; De Moortel et 
al. 2000) and Hinode (e.g. Ofman & Wang 2008) confirms 
the existence of magnetohydrodynamic (MHD) wave motions 
in the coronal plasma (see also reviews by De Moortel 2005; 
Nakariakov & Verwichte 2005; Ruderman & Erdelyi 2009; 
Goossens et al. 2011; De Moortel & Nakariakov 2012). MHD 
wave theory suggests that the corona potentially supports a wide 
variety of distinct classes of wave motions, namely the Alfven 
wave, and fast and slow magnetoacoustic waves. However, the 
applicability and appropriateness of such classifications has re- 
cently provoked intense discussion, for instance reports of hav- 
ing observed Alfven waves in the corona (Tomcyzyk et al. 2007; 
De Pontieu et al. 2007) are contested by, e.g., Erdelyi & Fedun 
(2007) and Van Doorsselaere et al. (2008). 

The 'classic' terminology of MHD waves originates in 
the analysis of modes supported by magnetically-unidirectional, 
homogeneous plasmas of infinite extent which consider either 
plane wave solutions or utilise the method of characteristics 
(viz. Riemann decomposition) (examples of such analyses can 
be found in, e.g. Friedrichs & Kranzer 1958; Lighthill 1960; 
Cowling 1976; Goedbloed & Poedts 2004). Here, three distinct 
modes are permitted, the Alfven mode, and the fast and slow 
magnetoacoustic, and their behaviour and nature is well under- 
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stood. For a low-jS plasma, it is found that the Alfven wave is 
a transverse, purely-magnetic wave, propagating at the Alfven 
speed, guided by the magnetic field. The fast magnetoacoustic 
wave is found to propagate roughly isotropically at the fast speed 

(cf = ^jc 2 A + c 2 s , Cf w ca where <c 1), and can travel along 

and across magnetic fieldlines. The slow magnetoacoustic wave 
propagates longitudinally along the magnetic fieldlines, roughly 
at the sound speed. 

However in the solar corona, effects including gravitational 
sttatification, inhomogeneous density profiles and multiple- 
source magnetic field geometries call into question whether 
these three, classical modes are still valid. The departure to inho- 
mogeneity typically introduces a variety of phenomena such as 
resonant absorption (in the corona see, e.g., Ionson 1978; 1982; 
1983; Hollweg 1984; Ruderman & Roberts 2002) and phase 
mixing (e.g. Heyvaerts & Priest 1983), all of which blur the dis- 
tinction between these separate modes (see also, e.g., Bogdan et 
al. 2003; McDougall & Hood 2007; Sousa & Cunha 2008; Cally 
& Hansen 201 1; Hansen & Cally 2012). 

Typically theorists preserve the concept of the three MHD 
modes and introduce the concept of conversion and coupling be- 
tween the constituent modes as waves encounter certain inho- 
mogeneous features. However, as highlighted by Goossens et al. 
(2011), we must take care to understand that realistic plasmas, 
snictly speaking, do not support three 'classic' wave modes 
but rather that a propagating MHD pulse encounters situations 
in which it assumes transient properties that qualitatively corre- 
spond to one (or more) of the homogeneous modes. Conversely, 
it should also be possible to locally determine hyperbolic char- 
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acteristics corresponding to the fast, slow and Alfven waves 
for any ideal MHD system via Riemann decomposition (for an 
overview see, e.g., Goedbloed et al. 2010) although, as discussed 
in Goedbloed & Poedts (2004), this is not always analytically 
possible for fully inhomogeneous cases (in particular for non- 
unidirectional magnetic fields). 

As such, it is unclear whether or not the classic MHD modes 
exist in the solar atmosphere as discrete entities or if such mode 
separation requires specific geometries. For example, as dis- 
cussed by Parker (1991), the existence of true Alfven waves as 
per Alfven (1942) is dependent upon special magnetic geome- 
tries which contain invariant directions. Despite this uncertainty, 
a theoretical framework which involves them conceptually is not 
necessarily obsolete. Mode interpretation and analysis that re- 
lates MHD wave behaviour in inhomogeneous, realistic plasmas 
to the comparatively simple, classic waves still has the potential 
to be a useful framework. 

In this paper, we are specifically interested with the be- 
haviour of MHD waves in the vicinity of the 3D magnetic null 
point (e.g. Parnell et al. 1996; Priest & Forbes 2000). Null points 
are topological features of coronal magnetic fields, predicted by 
magnetic field extrapolations such as Beveridge et al. (2002) 
and Brown & Priest (2001). At these magnetic null points, the 
magnetic induction is zero. Hence, approaching a magnetic null 
point, the coronal plasma is highly inhomogenous. In the so- 
lar atmosphere, null points have been identified as playing key 
roles in many processes, for example; in CMEs (in the magnetic 
breakout model, e.g. Antiochos 1998; Antiochos et al. 1999), 
magnetic reconnection (e.g. Priest & Forbes 2000) and oscil- 
latory reconnection (reconnection driven by wave-null interac- 
tions, see McLaughlin et al. 2009; McLaughlin et al. 2012) all 
of which are thought to play a role in coronal heating. As both 
waves and null points are ubiquitous in the corona (Close et al. 
2004, Longcope & Parnell 2009 and Regnier et al. 2008 give a 
rough estimates of 1 .0 - 4.0 x 10 4 null points) wave-null interac- 
tions are inevitable and thus are arguably a fundamental plasma 
process in the solar atmosphere. Hence, our present research is 
specifically concerned with the extension of classic MHD wave 
theory about null points, which in a broader sense contributes 
to both the theory of MHD wave behaviour in inhomogeneous 
media and the understanding of fundamental aspects of coronal 
physics. 

MHD wave behaviour in the neighbourhood of magnetic null 
points have been extensively studied in 2D models. Bulanov 
& Syrovatskii (1980) performed the first investigation of MHD 
behaviour about a 2D null and noted that in a 2D geometry 
the motions governing the Alfven mode and the fast magnetoa- 
coustic modes are decoupled, permitting analysis that considers 
the modes separately. The transient features of fast and Alfven 
waves in various 2D null point geometries within a/3 — plasma 
were studied extensively in a series of papers by McLaughlin 
& Hood (2004; 2005; 2006a). Again, these authors found the 
two modes decoupled and identified key propagation features for 
each mode. The fast wave propagates isotropically (at the same 
characteristic speed both across and along magnetic field lines), 
with behaviour dictated by the Alfven-speed profile, propagating 
from regions of high to low Alfven speed resulting in a refrac- 
tion effect which focuses the wave energy into the null point. 
Meanwhile, the Alfven wave is confined to follow magnetic 
fieldlines, thus leading the wave energy to accumulate along the 
separatrices. Due to the resulting current build-up in these re- 
gions, these papers concluded that magnetic null points are likely 
locations for localised heating events in the corona. f3 + and 
nonlinear behaviour has also been investigated by McLaughlin 



& Hood (2006b) and McLaughlin et al. (2009) respectively. A 
comprehensive overview of the whole topic is given in the re- 
view paper of McLaughlin et al. (201 la). 

Of course, 2D models only give an initial grounding in the 
physics of realistic null points and for a full understanding we 
must turn to 3D models (as singularities, null points are in- 
escapably 3D and any 2D X-point configuration in fact only cap- 
tures the physics of a null line of infinite extent). It is not clear 
to what extent the characteristics and behaviour of waves about 
2D null points transfer to the fully 3D case, and surprisingly 
few papers have been written that address MHD wave behaviour 
about a fully 3D null point. Most papers have focused on the 
dynamics of current accumulation over time (in an attempt to lo- 
cate regions where reconnection is most likely to occur) rather 
than focus on the transient propagation features of the individual 
modes. 

Galsgaard et al. (2003) consider a proper, /3 = null point 
and introduce a twist wave (what they call an Alfven wave) 
which is generated about the spine and eventually accumulates 
on the fan (perhaps analogous to Alfven waves behaviour around 
2D null points). In addition, they observe a small amount of cur- 
rent accumulation at the null itself, which they suggest is due to 
nonlinear generation of a fast wave. No linear coupling between 
wave modes is observed, and in their linear analysis they find 
that the wave equations for the fast and Alfven modes decouple. 
However, this is not surprising since their proper 3D null has az- 
imuthal symmetry, and so the system is actually 2.5D, not fully 
3D. 

Pontin & Galsgaard (2007), Pontin et al. (2007) and 
Galsgaard & Pontin (2011a; 2011b) performed numerical sim- 
ulations in which the spine and fan of a proper 3D null point 
are subjected to rotational and shear perturbations. They found 
that rotations of the fan plane lead to current density accumula- 
tion about the spine, and rotations about the spine lead to current 
sheets in the fan plane. In addition, shearing perturbations lead 
to 3D localised current sheets focused at the null point itself. 
Again, this is in good agreement with what we may expect for 
MHD wave behaviour from the 2D studies, i.e. current accumu- 
lation at specific, predictable parts of the magnetic topology. 

The first study of MHD wave behaviour in the neighbour- 
hood of a fully 3D null point was investigated by McLaughlin 
et al. (2008). These authors examine the fast and Alfven waves 
about both proper (i.e. 2.5D) and improper (fully 3D) null points. 
The authors utilise the WKB approximation to determine the 
transient properties of the modes in a linear, ft — plasma 
regime. Their findings strongly suggest that the features of MHD 
waves about 3D nulls are not that different to the 2D results: the 
fast wave propagates across magnetic fieldlines according to the 
Alfven-speed profile and the Alfven wave is confined to mag- 
netic field lines, leading to the waves accumulating in particular 
topological regions of the null point. However, their implemen- 
tation of the WKB method is unable to address the question of 
whether modes couple in the fully 3D geometry, since their first- 
order WKB solution implicitly precludes the possibility and con- 
strains the waves to see the magnetic field as locally uniform. 

Thus, as it stands the question of what is the true behaviour 
of MHD waves in the neighbourhood of fully 3D magnetic null 
points remains unanswered. Specifically, what is the nature of 
the propagation and evolution of each MHD wave and, critically, 
what is the efficiency of mode-coupling and/or conversion due to 
the magnetic geometry and/or due to nonlinear effects? 

This paper investigates the behaviour of fast magnetoacous- 
tic waves about fully 3D null points via numerical simulations 
and looks for evidence of coupling to the Alfven mode in any 
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Fig. 1. Left: The azimuthally symmetric proper 3D null point (e = 1). Right: An improper 3D null point where fieldlines are predominantly aligned 
with the x-axis (e = 0.5). 



form (e.g. linear, nonlinear and due to the 3D magnetic geome- 
try). The paper is structured as follows: §2 outlines the specifics 
of the model and is subdivided into sections on the structure 
of the null points considered (§2.1), the mathematical model 
(§2.2), the method for isolating individual wave modes (§2.3), 
the specifics of the numerical solution (§2.4) and a brief outline 
of the supporting WKB solution (§2.5). §3 and 4 present the re- 
sults of the simulations for the improper and proper radial null 
points respectively, and the conclusions are presented in §5. 



2. Mathematical Model 

2. 1 . Magnetic Null Points 

In this paper, we investigate MHD wave behaviour about 3D po- 
tential magnetic null points. In Cartesian coordinates, these null 
points take the form: 



B = [x,ey,-(e+l)z] 



(1) 



where the parameter e is related to the eccentricity of the mag- 
netic field lines, controlling the direction in which they predom- 
inantly align. Parnell et al. (1996) investigated and classified the 
different types of linear magnetic null points that can exist (our 
e parameter is called p in their work). For equation (1), the null 
point itself, i.e. B = 0, occurs at the origin. 

Priest & Titov (1996) identified two key topological features 
present near any 3D null point: the spine and the fan plane. The 
spine is an isolated fieldline that approaches, or leaves, the null 
point along the z-axis. In this paper, we only consider e > to 
restrict our attention to positive null points, and as such the spine 
represents fieldlines approaching the null from above and below 
the z — plane. Thus, the z — plane, known as the fan, consists 
of radial field lines, confined to the z — plane and pointing 
radially away from the null point. 

The eccentricity parameter e alters the field topology as fol- 
lows: 

• For e = 1 , the magnetic null point has azimuthal symmetry 
about its spine, with no prefered direction for fieldlines, and 
is known as a proper null. Thus, as mentioned above, this is 
actually a 2.5D null point. 



• Null points which deviate from this cylindrical symmetry are 
known as improper nulls. For < e < 1 the fieldlines curve 
to run parallel to the x-axis, and for e > 1 curve to run paral- 
lel to the y-axis. 

• For e = 0, we recover the simple 2D null point in the 
xz-plane, with a null line running through x — z = 0. 

See Parnell et al. (1996) for more comprehensive information on 
the classification of different types of 3D null. 

In this paper, we seek to determine the behaviour of fast 
MHD waves about potential 3D null points in the most general 
terms possible. However, analytical progress is unlikely when 
considering a completely general null point so we consider two 
different, specific null points that capture the range of topologies 
possible. The left panel of Figure 1 shows the magnetic field line 
structure of our first choice: the proper e = 1 null point (chosen 
for its azimuthally-symmetric fieldlines, which are (below) ex- 
ploited in analytical approaches). The right hand panel shows the 
second: an improper null defined by e = 0.5. Here fieldlines are 
predominantly aligned parallel to the x-axis, and by compari- 
son we address the question of how MHD wave behaviour alters 
with the field topology and eccentricity. 

2.2. Governing MHD Equations 

The three-dimensional, nonlinear, ideal, adiabatic MHD equa- 
tions are solved numerically for a /3 = plasma. In most of the 
solar corona the plasma-yS is very low, (J3 <K 1, Gary 2001) and so 
the f3 = approximation is arguably valid. Here, we are specif- 
ically concerned with addressing whether the 3D magnetic field 
about realistic null points introduces behaviour absent in pre- 
vious 2D models. Hence, in this study we utilise the f3 = ap- 
proximation intentionally to prohibit the introduction of the slow 
mode and restrict our attention to the behaviour of fast magne- 
toacoustic wave and its interplay with the Alfven mode, on the 
understanding that the results presented are a first step towards 
understanding how waves behave at realistic null points in the 
corona; which are both 3D and have ft + 0. 

Thus, the governing equations are as follows: 



<9v 

— + (v ■ V)v 

of 



VxB 



xB 



3 



dt 
dp 
dt 
dp 
It 
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V X (v X B) , 
-V-(pv) , 

-v-Vp — ypV • v . (2) 



Here, standard MHD notation applies: v is plasma velocity, 
p is plasma/thermal pressure, p is density, B is the magnetic 
field/induction, y — 5/3 is the adiabatic index, and p is the mag- 
netic permeability. 

We consider an equilibrium state of p = po, p = po (where 
po and po are constants), v = and equilibrium magnetic field 
B = Bo. Finite, small perturbations of amplitude a <K 1 are 
considered in the form p = po + oip\(r,t), p — po + ap\{r,t), 
x = + orv(r, t) and B = Bo + ab(r, t) and a subsequent nondi- 
mensionalisation using the substitution v = vv*,V = V*/L, 
B = -B B„, b = B b*, t = ft*, p\ = p p\ andpi = p p* is per- 
formed, with the additional choices v = Lfi and v = Bq/ yjppQ. 
The resulting nondimensionalised, governing equations of the 
perturbed system are: 



dv 

dt 

db 

dt 

dpi 
dt 

dp i 
dt 

N, 

N 2 

N } 
N 4 



(V x b) x B + Ni 

V x (v x B ) + N 2 
-V ■ v + N 3 

-y V ■ v + N 4 

dx 

(Vxb)xb-pi— -(▼• V)v 
dt 

V x (v x b) 
-V • (piv) 

v-(Vpi)-ypi(V-v) 



(3) 



where terms N, are the nonlinear components correct to 0(a 2 ). 
The star indices have been dropped, henceforth all equations are 
presented in a nondimensional form. The equations are merged 
into one governing PDE: 

^(Vx[Vx(vx Bo)]} x Bo + N 

N = {Vx[Vx(vxb)])xB„ 

+ {V x [V x (v x b)]} x b + (V x b) x [V x (v x B )] 
+ (V • v - v ■ V) (V x b) x B 
-pi{Vx[Vx(vxB )]}xBo 

- [(V x b) x Bo ■ V] v (4) 

The first term dominates the linear regime of O (a) and the terms 
N represent nonlinear terms of 0(a 2 ). 

This study utilises numerical solutions of the full nonlin- 
ear MHD equations (2), using the numerical method detailed in 
§2.4. In conjunction, analytical treatment of wave equation (4) 
derived via second-order perturbation theory is used to create a 
special coordinate system that will isolate and identify individ- 
ual wave modes (§2.3), create a WKB approximation (§2.5) and 
to seek analytical confirmation of observed phenomena in the 
numerical simulations where possible (§A). 
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Fig. 2. Fieldlines of magnetic induction Bo (blue) and the field C (red) 
perpendicular to both the magnetic induction and the flux function, 
in the case e = 1 for the y = plane. Due to azimuthal symmetry 
other such planes have identical fieldline structure. The flux function A 
(green) in this case is azimuthal about the spine, running antiparallel to 
above the fan plane and parallel to 6 below (® corresponds to a vector 
directed towards the point of view, and O away). These vector fields are 
everywhere perpendicular thus suitable for forming an orthogonal coor- 
dinate system, although fields C and A are undefined on the spine and 
fan plane. 



2.3. Isolating MHD modes 

This paper investigates the transient MHD mode behaviour 
about a general 3D null point as well as the possible mode cou- 
pling/conversion that may occur. To do so, it is necessary to 
construct a novel coordinate system that allows individual, pure 
MHD wave modes to be generated, identified and tracked. In 
turn, such distinction requires a clear definition of the possible 
modes that can occur. 

As outlined in § 1 , the concept and terminology of the three 
classes of MHD wave (Alfven, fast and slow) that originated 
within a homogenous MHD plasma description still remains 
valid and useful in describing inhomogenous media, although 
distinguishing between the individual modes can become in- 
creasingly difficult depending upon the exact nature of the inho- 
mogenity. The B — 0, uniform density, 3D magnetic null points 
under consideration in this paper are very specific examples of 
plasma inhomogenity, namely that of a nonuniform magnetic 
topology. Hence, the behaviour of the two linear waves sustain- 
able here (Alfven and fast) is solely due to the Lorentz force, and 
as such we find that the most appropriate and useful definitions 
of, and distinctions between, the two modes are purely in these 
terms. 

The nondimensionalised Lorentz force: 



F = (V x B) X B = (B ■ V)B - -(B • B) 



(5) 



is the sole driving force in our B — system, and is denoted F. 
The first term of the right hand side of the equation (5) can be in- 
terpreted as magnetic tension, and the second term is a gradient 
in magnetic pressure. We note that the force has no term parallel 
to the magnetic induction (F ■ B = 0) and it follows that no B = 
MHD wave can be driven by or associated with perturbations in 
this direction, i.e. they are all perpendicular perturbations rela- 
tive to total magnetic field. Thus, we define our two waves as 
follows: 
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1) The (linear) Alfven wave is driven only by magnetic tension. 
This results in a wave that propagates anisotropically as a 
transverse wave, confined to and guided along the magnetic 
fieldlines. 

2) The (linear) fast magnetoacoustic wave is driven both by 
magnetic tension and magnetic-pressure gradients. It prop- 
agates isotropically via a mixture of transverse and longitu- 
dinal motions (due to both magnetic tension and magnetic- 
pressure gradient respectively). 

These definitions are used for the specific purpose of creat- 
ing a special orthogonal coordinate system that can be used to 
distinguish between these two modes. Note that this will not 
be valid along the spine or along the fan plane, as here the 
waves are degenerate: at these locations both the fast wave and 
Alfven wave propagate only as transverse, tension-driven waves. 
Elsewhere, such a coordinate system can be obtained by exploit- 
ing a mathematically-invariant direction. 

For instance, if we consider initial perturbations to the equi- 
librium magnetic field (Bo) in an invariant direction denoted e,, 
then due to <9/<9e, = 0, the magnetic-pressure gradient is elim- 
inated. As such, this direction is inextricably linked with the 
Alfven mode as per our definition, i.e. perturbations in this di- 
rection give rise to waves driven only by magnetic tension that 
propagate as transverse waves along magnetic field lines. The 
existence of an invariant direction is just as much a physical re- 
quirement for the existence of Alfven waves as it is a convenient 
mathematical method for isolating them. Various authors have 
stated that in the absence of such an invariant coordinate, only 
the magnetoacoustic waves can exist, and hence only certain spe- 
cial topologies permit the Alfven mode (see, e.g., Parker 1991; 
Van Doorsselaere et al. 2008; Goossens et al. 201 1). 

Alternatively, we can consider perturbations in a direction 
that is perpendicular to both Bo and e, . Here both magnetic pres- 
sure and magnetic tension can be present and thus we assume 
this direction is associated with the fast mode as per our defini- 
tion. That is to say, both magnetic tension and magnetic-pressure 
gradients sustain an isotropic wave that propagates via longitu- 
dinal and transverse motions. 

To reiterate: in this paper we aim to construct a special co- 
ordinate system, for the 3D null, of the form Bo, e, and e, x Bo, 
such that motions in the e,-direction correspond to the Alfven 
mode and that motions in the e, x Bo-direction correspond to fast 
mode. 

Similar special coordiante systems have been previously 
utilised in various 2D MHD models, an example of which is the 
coordinate system used by McLaughlin & Hood (2004; 2005; 
2006a) to drive pure fast and Alfven modes at a 2D X-point. 
In fact, the essence of any 2D mathematical model is to impose 
an obvious invariant direction, and thus facilitating the creation 
of such a coordinate system. However, for a fully 3D null point 
(e + 0, e + 1), there is no obvious invariant direction and, unlike 
2D models, there is no clear way to discern which of the direc- 
tions perpendicular to Bo would be associated with unique MHD 
wave modes. 

We look for clues in the proper (e = 1) null case, which has 
obvious azimuthal symmetry. Considering the proper null point 
in cylindrical polar coordinates reduces the problem to 2.5D. We 
first seek to find the coordinate system which corresponds to 
unique wave modes in this specific case, then look for a gener- 
alised system (of which this would be a specific case) that would 
work with any 3D null point. 

The equilibrium magnetic field with e = 1 rephrased in 
cylindrical polars is Bq = [r, 0, -2z] where r 2 = x 2 + y 2 . Under 



our definition, the Alfven wave will be associated with the in- 
variant direction 0. Thus, it follows that the fast wave is asso- 
ciated with the direction that completes the orthogonal set, i.e. 
6 x Bo. Considering the linear components of the wave equa- 
tion (4) parallel to the orthogonal vectors 0, Bo and x Bo, in 
cylindrical polars for e = 1 leads to three decoupled wave equa- 
tions governing the linear Alfven, fast and slow (absent) mode 
respectively: 




This treatment offers a clear demonstration of (linearly) de- 
coupled waves for this specific case and provides confirmation 
that perturbations in will be associated with the Alfven wave 
and perturbations in x Bo with the fast wave. Such a treatment 
in cylindrical polars has been used in previous studies of proper 
magnetic null points, for example Galsgaard et al. (2003) utilised 
a symmetric helical motion (in our terminology, a perturbation 
in 0, to study an Alfven wave propagating about a proper, e = 1 
null point). However, it is not obvious how this method would 
transfer to the more general improper null case. 

2.3.1 . Special coordinate system: A, B and C = A x B 

The flux function/magnetic vector potential is defined such that 
Bo = V x A. For any simple, single potential null set-up, the 
magnetic helicity is zero (i.e. A • Bo = 0) and hence A is aligned 
with the invariant direction. Thus, we hypothesise that about po- 
tential null points the (linear) Alfven wave is always associated 
with the flux function direction. For the e = 1 null, we find the 
flux function, and thus the special coordinate system, would be 
(in cylindrical polars): 

Bo = [r,0, -2z] 

A = [0,-rz,0]oc£ (7) 
C = [2rz 2 ,0,zr 2 ]oc(£xB ) 

and we confirm that the vectors form an orthogonal set A • Bo = 
A • C = Bo • C = 0. Thus, to generate and isolate a pure Alfven 
mode, we consider perturbations in the A-direction only, and 
similarly for the pure fast wave, we must consider perturbations 
in the direction C, where C = A x Bo. The system for e = 1 can 
be considered a special case of cylindrical polars that changes 
about the fan plane, and is sufficiently similar to be consistent 
with the decoupled linearised wave equations for the cylindrical 
polars (6). Figure 2 illustrates the vector fields for the e = 1 case. 

By determining the flux function for general e, we extend 
this to create an orthogonal, curvilinear coordinate system that 
is based on the equilibrium magnetic field (Bo), the flux function 
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(A) and their cross product (C = A x Bo). These general vector 
fields are found to take the Cartesian form: 

B = [x,ey,-(e+l)z] 
A = [zy, -exz, (1 - e)xy] 

C y =y[(l-e)x 2 + (6+l)z 2 ] 
C z = ez(x 2 +y 2 ) 

with unit normals B = B/|B|, A = A/|A| and C = C/|C|. Note 
that on the line of the spine (x = y = 0) and the fan plane (z = 0), 
the identity B = V x A no longer holds, rendering the system 
locally invalid. As such, the spine or fan cannot be used to drive 
pure modes under this coordinate system. However, there is a 
good physical explanation for this, due to a degeneracy between 
the fast and Alfven modes on the fan plane and spine, i.e. a wave 
propagating along the spine or along one of the fieldlines on the 
fan plane is driven by tension only, as a transverse oscillation, 
and so along these specific/particular fieldlines the fast mode is 
degenerate with the Alfven mode. 

This special coordinate system (equation 8) is at the heart 
of the work presented in this paper. It is used in the following 
manner: Initially, it is employed to introduce a pure, linear, fast 
magnetoacoustic wave at a computational boundary by driving a 
velocity pulse in the direction C. Numerical experimentation (§3 
& 4) confirm that this special orthogonal coordinate system suc- 
cessfully isolates the fast mode. The coordinate system is then 
used to seek evidence of linear and/or nonlinear coupling to the 
Alfven mode, i.e. if we drive a pure fast wave (v^ + 0, = 0) 
and mode conversion occurs then under our coordinate system 
velocity components in the direction A will become manifest. 

2.4. Numerical Solution 

Nondimensionalised versions of the governing nonlinear MHD 
equations (3) are solved using the LARE3D numerical code (a 
Lagrangian-Eulerian remap scheme for MHD, see Arber et al. 
2001) with magnetic equilibria corresponding to both proper and 
improper 3D null points. Each scenario introduces a pure, planar 
fast magnetoacoustic wave pulse along the upper z-boundary, by 
driving the following sinusoidal profile: 

v e = a sin (2nt) , v A = , v 6q = for < t < 0.5 (9) 

Numerical tests determined that the amplitude a = 0.001 was 
sufficiently small relative to other parameters to generate a lin- 
ear wave and for the rest of this paper we set a = 0.001. The 
other boundary conditions are set as zero-gradient conditions. 
The simulations utilise a uniform numerical grid with domain 
-1 < x < 1, -1 < y < 1, -1 < z < 2 and 720 x 720 x 1080 grid 
points, giving an effective resolution of 5x « 6y « Sz ~ 1/360. 
The results presented focus into the region -1 < x < 1, 
— 1 <>' < 1,-1 <z< 1, i.e. a subset of the full numerical 
domain. 

Since C is undefined along the spine and the fan, driving 
generates a planar wave with a small hole about the spine (for 
the top-boundary driven simulations) or a gap along the fan (for 
the side-driven simulations). This is essential to avoid the mode 
degeneracy on the spine and/or fan-plane fieldlines. 

After the near-planar pulse is generated, a small trailing 
wake is also introduced due to dispersion (e.g. Uralov 2003). 
Such a wake does not impact upon the dynamics of the main/lead 
pulse. 



2.5. WKB Approximation 

Semi-analytical WKB approximations of fast waves about both 
proper and improper null points were previously considered by 
McLaughlin et al. (2008). Succinctly, the WKB method is an 
asymptotic series approximation for low-frequency wave forms 
and has previously proven useful in understanding null point 
models. The WKB solution for fast waves is presented in 3 over- 
laying our numerical results. More information regarding the 
WKB method (Wentzel 1926; Kramers 1926; Brillouin 1926) 
can be found in Bender & Orszag (1978). 

However, McLaughlin et al. (2008) only consider the first- 
order WKB approximation, and thus their WKB implementa- 
tion cannot account for the possibility of mode coupling, which 
is one of the chief concerns of this paper. As such, comparison to 
the WKB method serves a two-fold purpose: if the solutions dis- 
agree it indicate possible occurrences of mode coupling or, alter- 
natively, it indicates numerical inaccuracies. In addition, due to 
our special coordinate system being undefined on the spine and 
fan, the near-planar waves presented in the numerical solutions 
have the aforementioned gap around the spine. However, the 
mathematical nature of the WKB solution allows one to avoid 
the question of whether it is possible to drive a pure mode on the 
spine or fan and simply determines which path a fluid element 
would take if it were possible. As such, discrepancies between 
the WKB solution and numerical solutions on the spine or fan 
plane are entirely to be expected and indicative of neither cou- 
pling nor numerical error. 

3. Improper Null 

Let us first consider the improper radial null point, in the specific 
case of e = 0.5, where a velocity pulse (equation 9) is driven in 
V£ along the top computational boundary (i.e. z = 2). The resul- 
tant propagation is seen in figures 3 and 4, which figures show 
|v| in the planes defined by both x = and y = respectively. 

In both figures 3 and 4 we observe a refraction effect 
which focuses the wave towards the null point. This refrac- 
tion/wrapping effect has previously been identified as a key fea- 
ture of the behaviour of the fast magnetoacoustic mode about 
2D nulls (e.g. McLaughlin & Hood 2004) and is dictated by the 
Alfven-speed profile around the null point, namely that the wave 
propagates from regions of high to low Alfven speed. 

Here, the x — plane (Figure 3) shows the wave propa- 
gating in the region with the 'gentlest' Alfven-speed profile: 

c a~I {y 2 + 9^ 2 ) an d tne y - plane (Figure 4) shows the wave 
propagating in the region with the strongest Alfven-speed pro- 
file: c 2 = i (4x 2 + 9z 2 ). This difference in Alfven-speed profile 
directly corresponds to the difference of the rate of refraction 
in these figures, i.e. we observe the y = plane (Figure 4) to 
"wrap" around the null point more readily than that in the x = 
plane (Figure 3). 

The WKB solution for the leading edge of a planar fast wave 
is also shown in Figure 3 and 4 (in white). The agreement be- 
tween the leading edge of the wave front observed in |v| and 
the WKB solution near the null point is very good: the (nu- 
merical) wave experiences refraction along the same profile and 
at the same rate as prescribed by the analytical approximation. 
Towards the boundaries the numerical front increasingly lags be- 
hind the WKB front as time evolves: this not a physical attribute 
of the fast wave but rather a consequence of the boundary con- 
ditions. However, this does not affect the key behaviour close to 
the null where the agreement between the two solutions is good. 
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Fig. 3. Contour plots illustrating the propaga- 
tion of |v| about the improper 3D null point 
for the x = 0, yz-plane. The position of the 
lead wave front according to the WKB solu- 
tion shown in white, which is in good agree- 
ment with the numerical results (sufficiently far 
from the boundary). 



There is a further discrepancy where the WKB solution exists 
across the spine whereas the pulse does not. However, this is ex- 
pected as we purposefully did not drive any component along 
the spine due to the mode degeneracy discussed in §2.4. 

The agreement between the numerical solution and the WKB 
solution (which itself precludes mode coupling) suggests that 
the pulse is propagating only as the fast magnetoacoustic wave 
and thus no conversion to the Alfven mode occurs. To verify 
this, we now utilise the special coordinate system to inspect the 
individual constituent modes of the pulse. 



3.1. C : proxy for the fast mode 

We now compare |v| and |vg|; viz. the whole disturbance and the 
constituent fast mode. We find that the spatial extent, i.e. posi- 
tions of the leading, middle and trailing wavefronts, of |vg] is 
identical to that of |v| at any given time, and that the difference 
in magnitude between the two is 0(a 2 ). Hence, we observe its 
(linear) transient behaviour about a proper null point to be domi- 
nated by the refraction effect, focusing the pulse at the null point. 
As the difference between |v| and |vj,| is non-zero, it is clear that 
nonlinear effects, possibly 0(a 2 ) mode conversion, are present 
but the two are linearly equivalent. We now investigate the other 
orthogonal directions. 



3.2. A: proxy for the Alfven mode 

By considering |vjj, we find that there is no velocity perturbation 
in this direction at any location or time in the simulation. Hence, 
no figure is presented, since |v^J is identically zero for all time. 
As detailed in §2.3, this velocity component is in the direction of 
the Alfven wave, and as such the results indicate that a pure fast 



wave propagating about an improper 3D null does not generate 
an Alfven wave due to either linear or nonlinear coupling. 

3.3. B : proxy for the field-aligned motions 

Figures 5 and 6 show the velocity component in |VgJ (motions 
aligned with the equilibrium magnetic field) over the same time 
scale as Figures 3 and 4. We see that the wave shown is a non- 
linear disturbance of 0(a 2 ) that, bar a small trailing wake, is 
confined to within the same spatial locations as the linear fast 
wave. However, this nonlinear disturbance is of a different ge- 
ometry to that of the fast wave and splits into multiple peaks on 
the approach to the fan plane. These multiple peaks appear to be 
akin to the lobes seen in McLaughlin & Hood (2006b). During 
the evolution, we find that the disturbance of 0(a 2 ) is entirely 
sustained by and dependent upon the propagating fast pulse of 
O(a), which acts as a progenitor wave that creates this daughter 
disturbance instantaneously as it propagates. 

As the fast wave propagates about the proper null point, 
it causes small disturbances to the equilibrium (which can be 
thought of as both a force free magnetic field and a convec- 
tive/flow free fluid) as it passes. This effect is the subsequent 
field/fluid response (as a consequence of frozen flux the two 
are unsurprisingly connected) that settles to the initial config- 
uration. Due to this dependency on the propagating fast mode 
(it eventually returns to zero once the fast wave has passed an 
area, as opposed to being generated but then propagating inde- 
pendently) we have used the term disturbance when describing 
this effect as opposed to wave or mode. The effect responsible 
for the nonlinear disturbance is related (but not identical to) what 
has been called the ponderomotive effect in MHD literature (see 
e.g. Verwichte 1999; Verwichte et al. 1999; McLaughlin et al. 
201 lb). As this clearly has had no effect on the mode conversion 
process (which is absent) this effect is not further investigated. 
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Fig. 4. Contour plots illustrating the propaga- 
tion of |v| about the improper 3D null point 
for the y = 0, xz-plane. The position of the 
lead wave front according to the WKB solu- 
tion shown in white, which is in good agree- 
ment with the numerical results (sufficiently far 
from the boundary). 



Note, the numerical experiment was repeated with the driv- 
ing of a pure fast wave along one of the faces perpendicular to 
the fan plane, i.e. a 'side-driven' fast wave. The key results were 
found to be the same as those of the top-driven simulation, con- 
firming that the coupling and transient behaviour is independent 
of spine/plane geometry. 

4. Proper Null, e = 1 

The numerical experiments were also considered for the proper 
null point (e = 1). Here, the results where found to be consistent 
with, and analogous to, those in §3, namely that: 

1) The whole pulse and the fast wave component, |v| and |vg|, 
are spatially identical and linearly equivalent, with a small 
difference of 0(a 2 ). In both, we find a pulse that refracts 
along the Alfven-speed profile, in good agreement with 
WKB solution. 

2) There is no coupling, linear or nonlinear, to the Alfven mode 
(|Ts | is permanently zero throughout). 

3) There propagating fast wave sustains an instantaneous, non- 
linear disturbance in |vg |. 

Additionally, as the proper null point assumes azimuthal symme- 
try analytical confirmations of mode decoupling can be obtained 
(see appendix A). The results only differ in the rates of refrac- 
tion and time for accumulation at the null point, as the proper 
null has a different Alfven-speed profile. 

5. Conclusion 

We have studied both the transient properties of fast magnetoa- 
coustic waves and the nature of linear and nonlinear mode cou- 
pling in the vicinity of two potential 3D null points: the proper 



(e = 1) null, which can be treated as a 2.5D problem using cylin- 
drical polars, and the fully 3D, improper (e = 0.5) null point. 
Regardless of the eccentricity, e, of the null point studied, we 
find that: 

• An initially pure, linear fast wave exhibits no linear or non- 
linear coupling to the Alfven mode in the neighbourhood of 
a 3D null point. 

• Due to the absence of coupling, the propagation of the fast 
wave is entirely dictated by the Alfven-speed profile which 
about a 3D null point leads to a refraction effect, focusing all 
the wave energy at the null point itself. 

• The propagating fast wave generates and sustains an instan- 
taneous and dependent nonlinear field-aligned disturbance. 

One of the chief aims of the simulations was to look for ev- 
idence of MHD mode coupling about fully 3D (i.e. improper) 
null points. Unlike the proper symmetric 2.5D null, it was not 
clear from analysis of the MHD equations as to whether wave 
modes about improper 3D nulls are coupled due to the asym- 
metric magnetic topology, nor whether these coupling terms are 
linear or nonlinear, i.e. it was not clear what the overall general 
effect of departing from azimuthal symmetry would have on the 
wave dynamics. To address this we have generated a pure fast 
magnetoacoustic wave in the vicinity of 3D nulls and analysed 
its subsequent propagation. We find that this wave exhibits no 
mode coupling to the Alfven mode at either a linear or nonlin- 
ear scale about both the proper (e = 1) and improper null point 
(e = 0.5). Further experiments, namely e = 0.25 and e = 0.75 
(values e > 1 have analogues in the range < e < 1, and e = 
recovers the 2D null point), retain these results. Thus, we find 
that lack of mode conversion is a general feature of all potential 
null points, i.e. magnetic fieldline eccentricity does not facilitate 
fast to Alfven mode conversion. 
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Fig. 5. The nonlinear velocity component |v j 1 
for the x = cut of the improper null point: a 
small disturbance is created that mostly occu- 
pies the same geometry as |vg|. Each contour is 
scaled relevant to the amplitude at the specific 
times to enhance the contrast of the features. 



We additionally find that, in both numerical experiments, 
the fast-mode pulse generates and sustains a field-aligned dis- 
turbance in |vb q |. Our analysis shows that this field-aligned dis- 
turbance is not self-sustaining: it is an instantaneous daughter 
disturbance resulting from the progenitor fast wave. It is evident 
that the disturbance in |vb | does not act as a mode-conversion 
mechanism in the cases considered here, since no signal is ever 
generated in |va|, and the disturbance has no feedback effect 
upon the main wave. 

Due to the absence of mode coupling, the key features of 
fast wave behaviour in the neighbourhood of the e = 0.5, im- 
proper null points are fundamentally the same as the behaviour 
in 2D, ft = null point studies (e.g. McLaughlin & Hood 2004), 
despite the inhomogeneous, fully 3D magnetic field. The dif- 
ference between propagating fast waves about different poten- 
tial null points is due only to the differing Alfven-speed profiles. 
Thus, the propagation of the fast wave is found to be entirely dic- 
tated by such profiles, causing the refraction effect which, over 
time, focuses all of the wave energy at the null point. This is 
confirmed by our numerical simulations and supports the con- 
clusions drawn from a 3D WKB approach in McLaughlin et 
al. (2008). Note, in the numerical simulations at large time, the 
pulse is so close to the null point that the resolution will even- 
tually become inadequate. Nonetheless, in our simulations, over 
a finite, resolvable time the wave energy accumulates in a small 
spatial region around the null point, resulting in an exponential 
steepening of current-gradients, hence resistivity will eventually 
become non-negligible, resulting in ohmic heating of the local 
coronal plasma via resistive dissipation (a conclusion directly 
carried-over from 2D studies, see McLaughlin et al. 2011a). 
Thus, we conclude that 3D null points are locations for pref- 
erential heating by passing fast magnetoacoustic waves. 

Wave motions and null points are ubiquitous in the corona 
and so these fast wave-null interactions are likely inevitable. The 



large body of theory, extended by the results presented here, 
indicates that this will result in localised heating events. Thus, 
the key question is does this make a signifigant contribution to 
coronal heating? To answer this, clear observational evidence 
for MHD waves around coronal nulls is needed, a calculation 
of resultant heating, and a survey of the prevalence of such 
events. Such observations would require the detection of coro- 
nal nulls using co-temporal high spatial/temporal imaging and 
magnetograms to study oscillations in the vicinity of coronal 
null points. It is possible that the Atmospheric Imaging Assembly 
and Helioseismic and Magnetic Imager aboard SDO may permit 
such a study. 

We cannot, however, assume the results presented here hold 
in the opposite case, i.e. a similar scenario for an initially pure 
Alfven wave, follows suit. In fact, results of Galsgaard et al. 
(2003) report this may not be the case. Such an investigation 
will be the subject of a future paper. 

Finally, we highlight that sufficiently close to the null point, 
as magnetic induction drops off, there will be a region where 
the magnitudes of the sound speed and Alfven speed become 
comparable: identified by Bogdan et al. (2003) as the magnetic 
canopy, viz. the ft — 1 layer. This was investigated for a 2D 
null point by McLaughlin & Hood (2006b) and has the effect 
of introducing the slow magnetoacoustic wave into the system, 
which is coupled to the fast wave (but both are still decoupled 
from the Alfven wave). Thus the cold, ft — plasma assumption 
does not completely capture the physics of MHD waves about 
real null points, which must contain a ft = 1 layer near the null 
point. Thus, to understand wave dynamics about null points in 
the corona, it is necessary to extend the model presented here and 
see how these results are modified when the ft = assumption is 
removed. 
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Appendix A: Analytical confirmation of decoupling 
at the proper null point. 

In our simulations we find no evidence of any mode coupling 
from the fast to the Alfven mode. In the case of the proper null 
point, as it assumes azimuthal symmetry, we can confirm this an- 
alytically by taking the 0-component of equation (4) considered 
in standard cylindrical polars for e = 1 : 
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Here, the linear terms 0(a) are shown on the first line, and 
nonlinear terms of 0(a 2 ) follow. For the proper null 6 oc A, and 
so we consider which acts as a proxy for A. For an initial dis- 
turbance that is purely a fast mode disturbance, we take bg and 
Vg as initially zero, yielding d 2 v g /dt 2 = for all time. Hence, 
equation (A. 1) states that if bg = vg = initially then vg = 0, and 
hence = for all time. This provides analytical confirmation 
that about the proper null there is no fast mode to Alfven mode 
conversion at either a linear or nonlinear level. 
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